Interface fluctuations in disordered systems: 
Universality and non-Gaussian statistics 
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We employ a functional renormalization group to study interfaces in the presence of a pinning 
potential in d — 4 — e dimensions. In contrast to a previous approach [D.S. Fisher, Phys. Rev. Lett. 
56, 1964 (1986)] we use a soft-cutoff scheme. With the method developed here we confirm the value 
of the roughness exponent £ ~ 0.2083e in order e. Going beyond previous work, we demonstrate 
that this exponent is universal. In addition, we analyze the generation of higher cumulants in the 
disorder distribution and the role of temperature as a dangerously irrelevant variable. 

PACS numbers: 46.65.+g, ll.10.Gh, 75.10.Nr 
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I. INTRODUCTION 

Quenched disorder plays a crucial role in a huge va- 
riety of physical systems. One of the most prominent 
examples for such systems is a domain wall in Ising-like 
systems. Such interfaces can couple to various types of 
disorder such as a random pinning potential which can 
be provided by impurities in the system. For interfaces 
with internal dimensions d < 4, such disorder is known to 
radically modify the structure of the interface and also its 
dynamics. The structure of such systems typically is self- 
affine, resembling pure systems at criticality. In contrast 
to the statics, the dynamics is very different from a criti- 
cal one, since it lacks the characteristic scaling relations. 
It is dominated by high energy barriers which lead to 
an exponentially slow, glassy dynamics. This behavior 
was observed experimentally by position-space imaging 
techniques^ and it should also be possible to study such 
systems by means of scattering techniques .□ 

The theoretical understanding of such systems has 
made substantial progress during the last 15 years. Al- 
though the self-affinity of the structure suggests the use 
of renormalization-group (RG) techniques to calculate 
the characteristic exponents (such as the roughness ex- 
ponent), the analysis intricate since the flow of an infinite 
number of relevant parameters has to be studied. This 
can be achieved in the framework of a functional renor^ 
malization group which was employed by D.S. Fishem 
on the basis of the hard- cutoff renormalization group 
(HCRG) scheme of Wegner and Houghtonjj This study 
was the prototype for subsequent generalizations to a va- 
riety of different physical systems with a higher num- 
ber of components of the displacement fieldtl or to peri- 
odic systems&El. The latter class includes systems such as 
charge-density waves, Wigner crystals, and vortex lines 
in type-II superconductors. In all these systems disorder 
plays a qualitatively similar role. 

However, this HCRG is known to suffer from patholo- 
gies related to the sharp cutoff as noted already in Ref. 
||. This cutoff procedure leads to a long-ranged and os- 
cillatory behavior of the field correlations in real space, 



which leads to the generation of highly nonlocal inter- 
actions on a coarse-grained level. It requires particular 
care to interpret these interactions as the renormalization 
of local quantities such as the interface stiffness. These 
pathologies emerge not only in the context of disordered 
systems but also in pure systems, for example in the sine- 
Gordon model which describes the roughening transition 
of crystal surfaces. In the latter context, it was pointed 
out by Nozieres and GalletH that these pathologies may 
affect the renormalization of the stiffness constant and 
ultimately the scaling behavior at criticality. 

For this reason it is of principal interest to reinvestigate 
the prototype model for disordered systems in the frame- 
work of a soft-cutoff renormalization group (SCRG). In 
this scheme fluctuations are regularized on small length 
scales by a smooth cutoff function. There exist sev- 
eral variants of SCRGs (for recent review articles on RG 
schemes* see e.g. Refs. ^,^0|). The scheme of Wilson and 
KogutlB is a SCRG. Unfortunately, it is very clumsy to 
work with since already the free interface is described by 
a rather complicated-Gaussian fixed point. We use the 
scheme of Polchinskit3, which we find convenient for our 
purposes. 

Starting from this SCRG method we confirm the valuft 
of the roughness exponent found in the HCRG scheme.cl 
In addition, we explicitly demonstrate its universality, i.e. 
its independence from the cutoff function. In this respect 
our calculation parallels the demonstration of univetsaL 
ity for the bulk exponent rj of the pure Ising model.E3c3 

In most theoretical analyses disorder is assumed to be 
Gaussian distributed, i.e. that higher cumulants of the 
disorder distribution are negligible. Within the present 
scheme higher cumulants play a central role. Although in 
general these cumulants are non-universal, we determine 
their functional form at the RG fixed-point point since 
they are physically meaningful for energy fluctuations on 
large scales and since they will be needed for the analysis 
of the model to subleading order in e. 

In the next section || we specify the model. In order to 
keep this paper self-contained we include a brief scaling 
analysis in Sec. pill. The application of the Polchinski 
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SCRG scheme to the present model is discussed in Sec. 



IV 



The renormalization flow equations are derived and 
evaluated in Sec. |v| and our results are discussed in Sec. 
VI. 



II. MODEL 



Our analysis is based on the model Hamiltonian 



n 



J d d x{l(V^ + V(t( X ),x)}, (1) 



which is composed of an elastic energy and a pinning en- 
ergy. 7 is the elastic stiffness constant of the interface. 
V represents a quenched random potential acting on the 
interface. In the simplest case the disorder potential V is 
Gaussian distributed with zero average and short-ranged 
correlations in its dependence on x and <f>, 



y( ( /> 1 ,x 1 )T/(0 2 ,x 2 ) = Bfa - 2 )<5 (d) ( x i - x 2 ). (2) 

In principle one could allow for a more general form of 
the correlator (g) , where the dependence on x and <fi does 
not factorize and which has a finite correlation length also 
in x directions. However, in a coarse-grained description 
this correlation length shrinks to zero and one ends up 
with a dependence of the form (||) . Therefore this corre- 
lation length can be taken as zero right away without a 
modification of the large-scale properties of the system. 
On the contrary, the </> dependence of the correlator has 
to be described by the function R((f>) with a finite width. 
Otherwise the typical force density — J^V(4>, x) would be 
infinite and the problem would be ill defined. 

We subsequently focus on disorder of the random- 
potential type, where R{(j>) decays for large values of (f>. 
Thus, we exclude systems with random-field disorder, for 
which many large-scale features can he, derived already 
from Imry-Ma type scaling arguments. t3 We further spe- 
cialize to short ranged random potentials where R{4>) de- 
cays faster than any power of 4>- I n addition, we assume 
a statistical reflection symmetry <fi ~ * — <t> such that the 
correlation function is even, R(— 4>) = R(4>). 

For the analysis of the .model it is convenient to apply 
the standard replica trickED in order to anticipate the dis- 
order average and to restore translation symmetry. After 
replicating the system n times and averaging over disor- 
der, one obtains the Hamiltonian 



(3a) 
(3b) 



W — H c \ + Hpin, 

n „ 

n * = lY, ddx ( v< n 2 > 
z a =i j 

i n r 

«pin = -^ £ d d xR(r{x)-^(x-)), (3c) 

which we have decomposed into an "elastic" part and a 
"pinning" part. In this formulation, the quenched dis- 
order is represented by an effective interaction between 



different replicas. T denotes the temperature of the sys- 
tem. 

Since we have taken the correlation length in x di- 
rection as zero, the interaction is local in x and the in- 
teraction energy is invariant under an arbitrary replica- 
independent tilt Q (x) — > (j> a (x) + S(f>(x) of the interface. 
This symmetry will play an important role below. 

As usual, the partition function can be written as a 
functional integral 



Z = Tr d 



(4) 



with Tr = n k IL / #"(k) as integral over the Fourier 
modes <fi(k) = J d d x e~ lk ' x </>(x). This model has to be 
regularized at large momenta by a momentum cutoff A 
related to a microscopic length a = 2t:/A. Our particu- 
lar ch oice for the regularization will be discussed in Sec. 
IV B. 



III. SCALING ANALYSIS 

Before we study this model by means of the renor- 
malization group technique it is instructive to perform 
a scaling analysis. Although a lucid presentation of this 
analysis can be found in Ref. |^, we give a brief summary 
thereof in order to make this article self-contained. 

To describe the shape fluctuations of the interface we 
are primarily interested in the pair correlation 



(0(k)0(k')) = lim(0 a (k)^ tt (k')) =: C(k)<J(k + k') 

n—*0 



(5) 



(here a is not summed over) after averaging over thermal 
fluctuations (denoted by (• • • )) and over the disorder dis- 
tribution (denoted by 7TT ). Since the latter average has 
been anticipated in the replicated system (H) it no longer 
appears in the central term in Eq. (|E|). 

On large length scales k -C A the displacement corre- 
lation is found to be self-afflne, 



C(k) oc -k' d ~ 2 < 

7 



(6) 



with a roughness exponent £. For £ > 0, the interface is 
rough and the relative displacement 



w(x-x') := ([«£(x) 



x' 



= lim([f(x)-f(x')] 2 ) 



increases like 



n— >o 



u>(x) 



,2C 



(7) 



(8) 



on large scales x 3> a. We always assume £ < 1 since 
otherwise the model breaks down since higher powers of 
V</) become relevant in the elastic energy. For £ < 0, the 



2 



interface is flat and u>(x) converges to a finite value w c 
for x — > oo and 



w(x) 



,2C 



(9) 



In order to examine the relevance of disorder we an- 
alyze the properties of the model under a rescaling of 
lengths, field, and temperature: 



x = e xt, 
6 Q (x) = e« #(x<), 
T = e 6t To. 



(10a) 
(10b) 
(10c) 



Here we have introduced the scaling parameter t - which 
can be viewed as logarithmic length scale - and the en- 
ergy scaling exponent 9. The role of 0, which is ab- 
sent in usual critical phenomena, will become clear soon. 
The scaling hypothesis requires the statistical weights 
and therefore Ti/T to be invariant under rescaling. The 
elastic energy remains invariant only if the stiffness is 



rescaled according to 
differential form 



D (d-2+2(-9)£ 



7, which reads in 



d_ 

dt 



- /e = (d-2 + 2C- 0)7/. 



seal 



Thus the stiffness is scale invariant for 



9 = 2( + d- 2. 



(11) 



(12) 



a{3 



to 

AT 



r (14) 



which is the contribution to (||) bilinear in the displace- 
ment. 

If we reexamine the relevance of the couplings 
tom with the "random-force" exponents, we find 
■jl\ scal R2m.i = (4 - d)(m - l)R 2m ,i- Thus the couplings 
tom with m > 1 are irrelevant only in d > 4. Their rel- 
evance in d < 4 indicates that the roughness cannot be 
obtained from perturbation theory and that the correct 
value of the roughness exponent most likely is not given 

by Crf- 

In order to obtain the roughness exponent in d < 4, a 
renormalization group (RG) calculation is required. This 
RG has to be a functional RG since all terms in the Tay- 
lor series of R((j>) are relevant, as £ > and the relevance 
of coefficients i?2m increases with increasing m for any 
value of ( > 0. 

Before we start such a RG calculation, it is worthwhile 
to mention that a good estimate of the scaling expaoppt 
in d < 4 can be obtained from the Flory argument .Ejl3 
In this argument one assumes that for a rough inter- 
face the short-ranged disorder correlator can be approxi- 
mated by R(^>) ~ R*5(<f>) with a weight R* := J d<j> R{4>). 
Then ^| scal i?; = (d - C ~ 20) R^. The requirement that 
|| j< = and || ii| = then leads to the Flory 
value of the roughness exponent ( — Cf ■= ^jp- 



In the absence of disorder one can achieve the scale in- 
variance of both 7 and T with the choice 9 = and 
C - Cth — — ■ 

Since in the absence of disorder the interface is flat 
for d > 2 one may analyze the relevance of disorder 
(i.e. of 7ipi n ) performing a Taylor expansion of R((j>) = 
assuming analyticity of the (unrenor- 



V R 

Z-im>0 J2 



malized) correlator R (compare Ref. ||). A scale invari- 
ance of the statistical weights would require 



d_ 

di 



R 



2m, £ 



(d + 2m( - 29)R 



2m,e- 



(13) 



seal 



If we insert C = Cth an d 9 = 9 t h = into this equation, 
(d+m(2— d))R2 m ,t- In particular, 



Wefind di\ sc *\tom,l 

" ca i^2,£ = 2R 2 ,e and we conclude that disorder is rel- 

with m > 1 are less relevant 



evant. The couplings R 2rx 
for d > 2 under rescaling with the thermal exponents. 

The scale invariance of both TL C \/T and H P in/T can 
be achieved only for 9^0. Requiring ^| scal 7£ = and 

f e \ scal R2,i = one finds 9 = 0* =: 2 and C= Crf := 
which implies the roughness of the interface in d < 4. 
This "random-force" value of the roughness exponent 
is the value, one finds in a perturbative treatment of 
disorder , Ejc3 where the dependence of the pinning force 
on <j) is neglected an which is represented by the Hamil- 
tonian 



IV. RENORMALIZATION GROUP 

In d > 4 the couplings i?2 m are irrelevant for m > 1 
and the large-scale properties of the model are governed 
by the Gaussian random- force model ([l4|). In analogy 
to usual critical phenomena (see, e.g. Refs. |il"| , p3"| ) we as- 
sume that inci = 4— e<4 dimensions the large-scale fluc- 
tuations are described by a fixed point which is close to 
this Gaussian fixed point for small e. However, unlike for 
usual critical phenomena, this will b e a zero-temperature 
fixed pointQ According to Eq. (10c), which is equivalent 
to 



d_ 

di 



T,= 



(15) 



seal 



the effective scale-dependent temperature Tg flows to zero 
on large length scales x oc e e provided disorder increases 
the roughness of the interface, i.e. C > Cth and 9 > 
according to Eq. (|lT 



A. Rescaling 

Since we must perform a functional RG analysis, we 
first reformulate the flow of the couplings (such as 7 and 
R 2 m in sectieri HI) resulting from the rescaling (10) in a 
closed formBil 



3 



d_ 

d£ 



-0T' 



dTT' 



(x) 
(16) 



Here and henceforth we drop the subscript £ for simplic- 
ity of notation. In writing ^| scal we stress that this is 
only the scaling contribution to the flow. To obtain the 
full RG flow, this contribution has to be combined with 
a contribution ^|. which arises from integrating out 
modes of the field <b. 



C. Mode integration 

In Polchinski's schemejll the mode integration works 
as follows (see also Ref. ^5|) . One introduces an additional 
field d> < in the partition sum and defines d> > := <p — d> < . 
This can be achieved in such a way that 



xV e -*W > - G > 1 * > )-^p-[* < +* > ] (20) 



B. Regularization 

To calculate this second contribution, Fisher! and Ba- 
lents and Fishera have used the "hard cutoff" scheme 
of Wegner and Houghtoncl in their analysis of the 
present problem- We choose the "soft cutoff" scheme 
of PolchinskiE3Ej for the reasons outlined in the intro- 
duction. 

We regularize the theory by modifying the propagator 
using Schwinger's proper time method (see, e.g., Ref. 25). 
To this end, we rewrite 7i e i in Fourier space, 



H e \ = -(d>, G 



It, 



(17) 



Herein the propagator 

G A (k) := G oc (k)/(fc/A) with Goo(k) = l/( 7 k 2 ) (18) 

is regularized by the cutoff function /. In order to sup- 
press fluctuations on short length scales, f(z) has to van- 
ish for large z, whereas it has to satisfy / — > 1 for z — > 
in order not to modify the properties of the model on 
large length scales. f(k/A) can roughly be interpreted 
as "weight of modes" : In the calculation of the pair cor- 
relation (^) one may consider the density of modes in 
Fourier space as being reduced to a fraction / with re- 
spect to the unregularized model. 

We will keep the function / general as long as pos- 
sible, which is desirable to verify the independence of 
the results on the regularization procedure. However, we 
implicitly assume / to be monotonous in order to have 
smooth functions and to have no further intrinsic length 
scales. For illustrative purposes we occasionally choose 
the specific form 



-z 2 /2 



(19) 



As long as / _1 (fc/A) is analytic, a Taylor expansion of 
/ _1 (fc/A) shows that the regularization is achieved by 
contributions to the Hamiltonian which are irrelevant on 
large length scales since they involve higher orders of k. 
The HCRG method of Wegner and HoughtorD can be 
considered as non-analytic choice f(z) — 6(1 — z). 



The "slow modes" d> < are regularized by the propaga- 
tor G< := G\ K with an infinitesimally reduced cutoff 
A< := e~ A, whereas the "fast modes" d> > have a prop- 
agator G> := Ga — G<. In Eq. ( po| ) — which is derived 
from Eq. (Q) in Appendix ^ — proportionality factors 
independent of d> have been dropped. 

In the representation ( p0| ) the modes (f> > can be inte- 
grated out. This can be done exactly in the limit, where 
G> is small (it is of order d£) since then also d> > is small 
(of order (d£) 1 ^ 2 ). Then an expansion of 'H v \ n [4> < + d> > ] 
to second order in d> > is sufficient to establish the differ- 
ential RG. For infinitesimal d£ the propagator of the fast 
modes is 



G>(k) = g(k)de + 0{d£ 2 ), 
9(k)=AA GA (k). 



(21a) 
(21b) 



Integration over the field d >> then yields a flow (see Ap- 
pendix |a|) : 



d 
~d( 



512 



T 



(22) 



We introduce the abbreviations f. ., := f. f. f. . . . , 

J t := J d d Xi, 4>i := 4>(Xi) and := gfa ~ Xj) to keep 
expressions compact. We call the first term in Eq. ( |22j ) 
the "contraction" part of the generator (since legs of a 
vertex are contracted in a diagrammatic representation) 
and the second term the "composition" part (since new 
vertices are composed by linking two vertices). 

The flow contribution (|22| ) has a simpler form than the 
corresponding contribution in the HCRG scheme, since 
( p2| ) contains terms only up to first order in g (and second 
order in the Hamiltonian) , whereas the HCRG generator 
contains terms of arbitrarily high order. This difference is 
due to the fact that G> is a bounded function of order d£, 
in contrast to the HCRG scheme, where G> is singular for 
momenta at the cutoff, for which reason the power count- 
ing of orders in d£ breaks down and also higher powers in 
G> contribute to 0(d£). The simplicity of the generator 
is one important argument in favor of the SCRG scheme. 

The complete functional RG consists of the combina- 
tion of integration over modes (p2|) with rescaling (O) , 
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-H= - 

dl di 



H 



seal 



d_ 

di 



n. 



(23) 



Although we have defined the regularization in Fourier 
space, we ultimately find it more convenient to evaluate 
the RG flow in position space because of the functional 
structure of ^| int ^- 

To close this exposition of the RG method we empha- 
size that the RG (G3) is exact, provided the dependence 
of the functional Ti p \n[<t> < + 4> > ] on small <j> > (which is of 
order (d£) 1 / 2 because of Eq. (21a)) is captured in order 
di by a functional Taylor expansion to second order in 
c/) > . This analyticity assumption is common to HCRG 
and SCRG approaches and will be reexamined later on 
after we have found the fixed point for the disordered 
interface and the generation of nonanalytic features. 



D. General features 

We now turn to the evaluation of the RG flow for our 
model (||). Due to the locality of the disorder correla- 
tor in x space, the model owns a stochastic symmetry 
undcrj-ajtilt <fi a (x) — > cj) a (x) + v • x with a constant vec- 
tor v.ESo As a consequence of this symmetry there is no 
renormalization of the stiffness 7 arising from the mode 
integration. Hence the flow of 7 arises only from rescal- 
ing, 



d( 



7 = (d - 2 + 2C - 0)7. 



(24) 



This implies the validity of (|12|) not only in a scaling 
analysis but also at the RG fixed point. 

In contrast to 7, the couplings of the disorder part of 
the Hamiltonian, Ti. p - m , will not only be rescaled but also 
renormalized by the mode integration. This renormaliza- 
tion results not only in a modification of the second-order 
cumulant R of the disorder distribution but also in the 
generation of contributions to 7i p i n of an extended func- 
tional form. The functional form of the emerging terms 
can be recognized from the action of the generator (E^) 



on 7ipi n - The unrenormalized functional ( J3c| ) evaluates 
the field in two different replicas but only at a single point 
x. Therefore we call functionals of this type "2-replica" 
functionals and "1-point" functionals. The action of the 
first "contraction" term in ( p2[) results again in a 1-point 
functional. However, the "composition" term evaluates 
the field at two different positions, it is a 2-point and 
3-replica term. Fortunately, at the expense of the more 
complicated functional form we get "smaller" terms (i.e. 
of higher order in R). This 2-point term is of second 
order in R. Successive iterations generate p-point and 
(p + l)-replica terms in order R p . 

The extended functional form of the Hamiltonian can 
be captured in the form 



1 

T' 



2T 2 + 3!T 3 



where Q, 1Z, and S are 1-, 2-, and 3-replica terms, 
which represent cumulants of the disorder distribution. 
To become specific, let us denote the pinning energy 
of the replicated system before disorder averaging by 
TLv = Yl a Jx ^(^ a ( x )' x )- Then the contributions to the 
extended Hamiltonian can be identified as cumulants of 
TLv'- 



Q = (H V ), TZ = (H 2 V ) C , S = (H 3 V ) C ,. 



(26) 



(25) 



We will not keep track of Q since this is a field- 
independent constant because of the stochastic symmetry 
</>(x) — > 0(x) + 6<f> for an arbitrary constant 8<p. 

Anticipating that R — 0(e) at the fixed point, it is 
sufficient to retain p-point and (p + l)-rcplica terms to 
study the RG flow in 0(e p ). As we will see in Sec. [v], it is 
possible to retain the full functional form of the Hamilto- 
nian without need for truncations. To keep track of this 
functional form it is more convenient to perform the RG 
analysis in position space than in momentum space. 



V. RG FLOW AND FIXED POINT 

In this section we analyze the RG flow in order e 2 , from 
which we obtain the roughness exponent £ in order-ft- 
We find agreement with previous HCRG calculationscrtl 
in this exponent. In addition, we are able to demon- 
strate the universality of this exponent in the sense of 
its independence on the cutoff function /. We also ex- 
plicitly obtain the third-order cumulant S of the disorder 
distribution. n n 

In the HCRG analysis&EI it was not necessary to keep 
track of 3-replica terms for the determination of C in order 
e. Although 3-replica terms are generated, they do not 
feed back into the equation that determines £ (in order 
e) and there is no need to keep track of S. The situation 
is different in the SCRG scheme, since 1Z is renormalized 
only via S. However, 3-replica terms are generated also 
in the HCRG scheme. For the analysis of the consistency 
of the e-expansion one has to include also this term into 
consideration. 

In order to determine £ in order e, the it is sufficient 
to sufficient to consider the Hamiltonian in a functional 
form which can be parameterized by functions i?+ , 
and S according to 

+ W R ++ (cbf,<j>f;x 12 ) + -- - , (27a) 

a/3 Jl2 

S = J2 I S—itf ,C y ;xi2) + - , (27b) 

a J 12 

where we rewrite R+(<fi) = R(4>) to display that R is 
an even function of <f>. To compress notation we further 
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introduce 0"^ := 0"(xj) — 0^(xj). We thus keep func- 
tional contributions up to 3-replica and 2-point type. It 
is not necessary to keep these functionals in the most gen- 
eral form. Instead, the functions obey certain symmetry 
properties. R++ is an even function of both field argu- 
ments, S-- is an odd function of both field arguments 
(the detailed symmetry properties are specified below). 
Both functions evaluate the fields at two points and can 
therefore depend also on the distance xy := Xj — x 3 be- 
tween these points. 

A priori there is an ambiguity in representing a func- 
tional (such as TV) in terms of p-point kernels (such 
as the 1-point and 2-point functions R + and R++). 
For example, the simultaneous replacement i? + (0"^) — > 
bf) + Ai? + (0f ) and R+^f , J>f 



;xi 2 ) 



R ++ (4\ 13 , 02 /3 ; x i2) — Ai? + (0"^)5(xi 2 ) leaves the func- 
tional TZ unchanged. In this way one could absorb the 
1-point kernel into the 2-point kernel (in general, low- 
point kernels into higher-point kernels). In the case of 
S we avoid introducing a 1-point kernel. In the case 
of TZ we extract the 1-point kernel in order to ensure 
that the 2-point kernel contributes to TZ only in order e 2 
(see below). A unique distinction between the 1-point 
and higher-point kernels in a functional is achieved by 
requirement that the integral over higher-point kernels 
must vanish for a spatially constant field Q/3 (x) = 0q^ , 
in particular J d d x R ++ (4>^ , 0q^; x) = 0. 

The flow equation ( p3[ ) for the Hamiltonian can be rep- 
resented as flow equation of the parameter functions (for 
a diagrammatic representation, see Figjl]) 



df 



R, 



dt~ 



^R+{<t>) = {d - 20 + C0d}i?+(</>) + | J dtdaS—ttuO; x u ) ffu 
~ J dickS — (0i, </>i 5x12)312 H , 

^25x12) = {2d -26 + Qfadi + x, ■ Vi}R ++ (<pi, 2 ; x 12 ) - -did 2 S — (0i, 025x12)512 
+ ^S 12 J did 2 S — (01, 0i 5x13)513 H , 

ii, 2 ; xia) = {2d -39 + C0i5i + x, • V i}S__(0i, 025 x i2 ) - 3dfl+(0i)di? + (0 2 )5i2 + • • 



(28a) 

(28b) 
(28c) 



Vj := g^- is a partial derivative that acts only on the 
explicit position arguments xy but not on the arguments 
of the fields. In Xj ■ Vj summation over i is assumed im- 
plicitly (1 < i < p in p-point kernels). Further on, di 
denotes the partial derivative acting on the the ith field 
argument of the kernel function (which is not necessarily 
0i; if there is only one argument we drop the subscript, 
d := dr). 

In the flow equations ( |2S| ) we neglect all terms propor- 
tional to temperature since T flows to zero according to 
Eq. (|l5| ) . One immediately recognizes that even if there 
is only a function R + in the unrenormalized model, the 

kernel S is generated from the composition of two R+ . 

The kernel 5 then feeds back into the flow equations 

for i?+ and generates R++. 

Since R+(<fi) and g(x) are even functions, S-- has the 
symmetry properties 

S — (0i,0 2 ;x) = S — (0i, 025 -x) = S — (0 2 ,0i;x) 

= -5__(-0i,0 2 ;x) = -S , __(0i,-0 2 ;x). (29) 

These symmetries immediately imply 

R++ (01 , 02 ; X) = R++ (01 , 025 -X) = R ++ (02 , 015 x) 

= R++(-<f>i,<fa;x) = R++{<l>i,-<fa;x)- (30) 



The last terms in Eq. (28a) and Eq. (28b) are an 
"insertion of zero" since their contributions to TZ cancel 



exactly. They represent the shift of a 1-point term from 
the flow of R++ (0i, 02;xi2) to the flow of i£+(0i) in or- 
der to achieve J i? ++ (0 o , 0o; x) = according to the 
requirement specified above. This requirement has to be 

imposed on (unlike S ) in order to ensure that an 

e expansion can be performed consistently. 

In order to obtain the roughness exponent to lowest 
order in e, 



(1) 



h.o. 



(31) 



("h.o." stands for higher orders in e) we must determine 
the fixed-point Hamiltonian up to order e 2 , 



Wpin = eWg + e 2 wgl + h.o. 



(32) 



To lowest order there is only a 2-replica 1-point contri- 
bution 



n 



w = 

pin 2^ 



(33) 



In order e 2 the RG flow generates 2-point terms which 
are of 2-replica or 3-replica type: 
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n (2) = __ 



"pin 



-y 

a/3 

— y 
—y 

3!T 2 ^ 

a/37 



R 



i ) 



R 



12 



V 2 i( 



j(2) 



A Q 7. 



Xl 2 )- 



The solution of this differential equation (which we re- 
quire to preserve the rotation symmetry in x of the un- 
renormalized model) involves one constant of integration, 
which can be fixed by the requirement that the Fourier 
transform cr(k) be analytic at k = 0. This is a common 
requirement, which is crucial also for the solution of the 
(34) fa model, see Refs. p[Jll[. The Fourier transform of Eq. 
@ is (2 - e + k • Vk)cr(k) = g(k) which is solved by 



Inserting the Hamiltonian (^) and (|3J) into equation 
(^3|) we obtain the fixed-point equations for the parame- 
ter functions. We establish these equations and find their 
solutions in the order as they are generated by the RG. 

reads 



,(2) 



The fixed-point condition -^S__ '.< >;, 

= {2 + x. i - V 4 }si 2 l 

-3dR^\fa)dR^\fa)g 12 . 



2;xi 2 j 
i,02;xi 2 ) 



(35) 



(2) 

It is solved by a function S__((j)i, fa; X12) in which the 
dependences on all arguments factorize: 



S 



(2) 



92;Xi 2 ) 



-3dR<+\fa)dR<+\fa)<7(x 12 ). (36) 



The explicit spatial dependence is carried by the function 
a which is determined by the differential equation 



2cr(x) + x • Vcr(x) = -p(x) 



(37) 



a(k) = f 
Jo 



dz 



g(zk). 



(38) 



This solution is analytic at k = for d > 2 since g(k) 
is analytic. The explicit form of various appearing ker- 
nel functions is given in appendix [b| for the special cutoff 
function (|l9|). Since R^ is of order e, a is needed only 

in order e° to determine S^l_ from Eq. (|36"|). In gen- 
eral, g(k = 0) ^ and then also a(k = 0) ^ 0. The 

(2) 

k = contribution to the functional S__ could in prin- 
ciple be split off as a 1-point term. We refrain from doing 
so because both the 1-point term and the 2-point term 
are of order e 2 and this would unnecessarily increase the 
number of terms. 

In order e 2 the fixed-point conditions for R^' and R++ 
read 



= {1 -4( (1) + ( {1) fad}R { + ] (fa) + ^3 1 d 2 Si 2 !(<£ 1 ,0;x 12 )« ?12 -i Jd 1 a 2 S^.(<l>i 1 <h^vi)9v2, (39a) 
= {4 + x 4 • V,}7eg(^i,0a;xi 2 ) - -Si d 2 S^_(fa, fa; xi 2 )ffi2 + ^12 jf fcfc£i 2 l(&, fa; x 13 ) 5 i 3 - (39b) 



(39t), is found in the form 



Plugging the solution (B6J) into the fixed-point equation 



R ( +Ufa,fa;xi2) = -d 2 R ( +\fa)d 2 R { +\fa) P (x 12 ) (40) 



with a nonlocal kernel p satisfying the equation 

4p(x) + x ■ V j o(x) = c(x) - Cq5(x), 
where we defined 

c(x) := cr(x) 5 (x), 

co := / c(x) = c(k = 0). 

J X 



(41) 

(42a) 
(42b) 



The differential equation for p is solved in Fourier space 
in analogy to Eq. ( |37|) imposing the analyticity require- 
ment at small wave vectors, 



p(k) = r 

Jo 



co 



:(*k)]. 



(43) 



Since <7(k) is analytic, cr(k) and c(k) are analytic and so 
is p(k). It is important to notice that it is possible to 



find a solution which is analytic and finite for e — > only 



due to the presence of the last term in equation ( 39b). In 
the absence of this term, Co would be absent in Eq. ( |43| ) 
and pik = 0) cx 1/e would become singular for d — » 4^. 
This would contradict our assumption that the 2-point 
kernel contributes to 7i only in order e 2 ; instead it 
would contribute to 7i in order e and generate even more 
complicated contributions to H. in order e 2 . 

Now we turn to the determination of the fixed-point 

function Rv' from equation (39a). Inserting the given 

(2) 

solution for S__ into this equation we find 

= {1-4C (1) +C (1 Ud}R+ ) W 

+c [d 2 R { + ] {fad 2 R ( l ] {fa - 2d 2 R ( l ] (fad 2 R { + } (0)] . (44) 

Remarkably, the entire information about the cutoff func- 
tion / (as well as the surface stiffness 7) is contained in 
the constant Co- This constant can be eliminated by a 
rescaling R — > R/cq, after which our flow equation fipjn- 
cides with the one derived with the HCRG scheme.HH 
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Since the fixed-point equation Q44| ) is homogeneous, i.e. 
invariant under a simultaneous rescaling 

:= 60 and R+(4>) -> R+(4>) := fr 4 i?+(<W (45) 

for any real 6, it has a continuous set of solutions. 
Therefore one can search a solution of Eq. ( ]44| ) (which 
is a nonlinear differential equation of second order for 
R+^) without loss of generality with the initial condition 
R^iO) = 1. The second condition <9i?i 1} (0) = fol- 



The most distinctive feature of the fixed-point solu- 



lows from the symmetry of The fixed point has to 
be determined numerically. For a trial value of £W one 
integrates Eq. (|4|) from <j> = to <f> = oo. The fixed 
point describing elastic interfaces in disorder with short- 
ranged correlations is determined from the condition that 
R+\c/)) decays monotonously and faster than any power. 
We find 



C (1) = 0.208298(1) 



(46) 



and the fixed-point function and its derivatives are illus- 
trated in Fig. H As discussed in Ref. ||, from this flow 
equation (^3|) one actually can determine not only the 
fixed point and scaling behavior for short-ranged corre- 
lations but also for long-ranged correlations which we will 
not further analyze here. 



VI. RESULTS AND DISCUSSION 



The RG result (|3lJ) with (|6|) is close to the Flory 
value (f = f • For d = 1 the value £(d = 1) = | of 
the roughness exponent is known_jexactly.EZl In d = 2 
and d = 3 numerical calculations^ resulted in the val- 
ues ((d = 2) = 0.41(1) and ((d = 3) = 0.22(1). Al- 
though it is questionable whether the corresponding val- 
ues of e = 3 are small enough to justify the neglection of 
higher order terms, the lowest order estimate ( |3l| ) yields 
values ((d = 1) « 0.6249, ((d = 2) » 0.4166, and 
C(d = 3) 1=3 0.2083, which are surprisingly close to the 
exact or numerical values. 

Since / enters the fixed-point equation only via Cq 
which can be eliminated by rescaling, we have explic- 
itly shown universality, i.e. that Q 1 ' is independent of 
the cutoff function /. However, the fixed-point Hamil- 
tonian does depend on this function. In particular, also 
the disorder cumulants are nonuniversal. We wish to 
stress again that in our analysis the nonlocality has beea 
treated exactly, whereas FishertJ and Balents and FishertJ 
have truncated the flow to spatially constant fields which 
is sufficient to determine ■ The nonlocality is impor- 
tant for the evaluation of the RG flow in higher orders of 
e.EJE] However, only 5^1 feeds back into the flow equa- 
tion for R+ from which C' 1 ) is determined. In contrast, 

i?Y+ with the proper definition of p [where p(k = 0) = 
ensures that R ++ is of order e 2 ] does not feed back. Nev- 

(2) 

ertheless, it is important to kapw in order to deter- 
mine C in higher orders of e£3tHJ 



tion i*, the nonanalyticity of 



(i) 



as pointed out by 



Fisheru and Balents and FisherEI and which is implied 
by the fixed-point equation flli]) as follows. According 
to Eq. (||) a physically meaningful solution must sat- 
isfy R ( + } (0) > and R ( +\<t>) can be expected to de- 
cay monotonously decaying for increasing <j>. Search- 
ing a solution which is as smooth as possible, we may 
assume dR^'((f>) to be continuous and to vanish for 
(j) = 0. Evaluating Eq. at <j) — one then finds 

d 2 R { +\0) = -[(1 - 4C (1) )i?^ 1) (0)/c ] 1/2 which can have 
a solution only for £W < j. Then the second derivative 
of Eq. ( 44 ) implies 



9 3 < 1) (0 ± ) = ± V /(l-2C( 1 ))[-9 2 i?V 1) (0)/2c ]. (47) 

The sign in front of the root is determined by the re- 
quirement that the amplitude of d 2 R+\</>) should have 
a maximum at <p = 0. Thus d 3 R ( ^\o ± ) ^ 0, i.e. R^((j>) 

is discontinuous at <f> = 0. Consequently, 9 4 i?^(0) con- 
tains a singular contribution proportional to 5(<fi). 

Now we look back to verify that this nonanalytic- 
ity does not invalidate the RG analysis performed so 
far. Since at the fixed point S__ is proportional to 

dR+\<f>), discontinuities appear if S^]_ is derived twice 
with respe ct to one field argument. In the derivation of 

Eq. (28a) some terms [for example J x d 2 S (</>i,0;x)] 

have been dropped because of the symmetry properties 
( p9| ) of S This amounts to the implicit assump- 
tion d 3 R+(0) — 0, which seems arbitrary because of 
the discontinuity of <9 3 i?+. However, in general there 
are additional odd factors since the energy functional 
has to be even in <f>. In the aforementioned example 
f x <9!si 2 !(<£i,0;x) oc dR^i^dZR^iO), i.e. there is an 

additional odd factor dR+\cj)i). The corresponding con- 
tribution to the functional ( p4[) then vanishes after sum- 
ming over the replica indices irrespective of the value of 
d 3 R+(0). Therefore there is no need to retain this term. 

The situation becomes more severe as soon as a fourth 
derivative <9 4 i?+(0) appears. Since the fixed point func- 



tion R+\. involves second derivatives of R+' , this can 

(2) 

happen where second derivatives of R + + with respect to 
one field argument appear. In the above flow equations 
(p8|) and the resulting fixed-point equations in order e 2 
this is not the case. However, in these equations we have 
neglected terms proportional to temperature T. In par- 
ticular, we have neglected contributions to the flow of 

i,0 2 ;xi 2 ) cx 



(i) 



R++) which are proportional to T<9 2 i?' 2 



Td^R^t {cf>). From a similar contribution to the flow 

equation of d 2 R^ one expects d i R < ^\o) cx T _1 , where 
T is the effective temperature that decreases on large 
scales. Thus these particular terms are proportional to 
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T°. This means that temperature is a dangerously irrel- 
evant variable. (For a discussion of the role of tempera- 
ture in d = 1 see Refs. [3l]|}2].) However, since does 

not feed back into the fixed-point equation for Rv , the 
roughness exponent is not affected in leading order. But 
these terms give contributions to the fluctuations of free 
energy (that scale for large system sizes L proportional 
to L ) which do not vanish in the limit T — > + . 

To sum up, we have shown universality of the rough- 
ness exponent to leading order in e with the SCRG 
method elaborated here. Avoiding locality truncations, 
we have obtained the nonlocal functional form of the 
fixed-point Hamiltonian. In particular, we have deter- 
mined higher cumulants of the effective pinning energy 
distribution on large scales the implications of which will 
be examined elsewhere. cB 

Although in d = 2 and d = 3 the agreement between 
the RG result for £ and numerical values is quite good, 
it is of fundamental interest to examine the possibility 
to extend the analytic theory beyond this leading .-order. 
It was argued by Fishera and Balents and Fishera that 
because of the nonanalyticity of the fixed point function 
R+(4>) in 0(e) the next higher order after 0(e 2 ) to the 
fixed-point equations should be C(e 5 / 2 ). Consequently, 
the subleading contribution to £ would be of order e 3 / 2 . 
However, their reasoning is based on a T = argument. 
As shown above, temperature is not a truly irrelevant 
variable and the validity of a T = argument for T > 
is questionable. From preliminary studiesEj'EII we expect 
that the subleading contribution to £ is of order e 2 . The 
SCRG scheme presented here lays the foundation for an 
extended analysis of the functional RG beyond leading 
order in a systematic way. In addition, this method is 
free of the pathologies hampering the HCRG scheme and 
better tractable than other SCRG schemes. 
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APPENDIX A: MODE INTEGRATION 

In this appendix we sketch the derivation of Eq. (|2C 
from Eq. (EL and deduce the RG generator ( p2| ) following 
PolchinskiH 

Starting from G\ = G< + G> , elementary manipula- 
tions lead to the identity 



(Al) 



We consider the partition sum of the pure interface as 
functional of the propagator G, 



Z e \[G] := Tr^ e 



From Ga = G< + G> one can immediately derive 
Z c i[G<]Z c i[G>] 



-Z c i[Ga] = 



^[(G^ + G^ 1 ) 



-i, • 



(A2) 



(A3) 



Now the partition sum (|f|) is transformed by the fol- 
lowing steps: we use Eq. ([Al]), introduce the additional 
field </> < , regroup fields in the bilinear Hamiltonian, use 
([A3]) , and substitute 4> = < + 4 I> : 



2 = Xr^ e~ir^ ,G <^ < £) - T w pi»[0] 



1 



Tr Tr e -^(0 < ,(G- 1 +G > - 1 )0<) + i(0<,G-V)-^(0,G > - 1 0)-i-H pin [0] 



-2«i [Ga] 



Tr <Tr e~ & ^ ^ G <^)-M(4>-^)<G> L {^-^))-^H vi , 
Tr fA <Tr W) > e -2V(0 < ,G< 1 < )~2T(0 > ,G- 1 > )-iH pin [0<+0= 



-Z : ei[G<]Z e i[G>] 

2 c i[Ga] „ _!_,-.. , _!_, . ,.: 1 , i;, 

Z cl [G<]Z i[G>] 



2 c \[Ga] 

Mg<] 



Tr < e -5T w< ' G < 1<#,<) (e~T«Pin[0 < +^ > ]) 



(A4) 



M) apart from the field independent factors 



This is Eq. 
2 C \. 

The RG generator (^2|) follows from evaluati ng t he ex- 
pectation value in the last expression of Eq. ( |A4[ ) . The 
coarse-grained pinning Hamiltonian H' in [0 < ] for the slow 



modes is defined by 

and is found from a cumulant expansion: 



(A5) 
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= (^pin[0 < 4 


0> 




IT ^^P^ n 


0< 


A Taylor expansion of 7i p i n 
<P(xi) etc.) 




(H pin [0 < +< 


P > }}> =«pin[ 


0<] 




T 


(2 <5 


(H 2 pin l<l><+< 


= r / 

7 12 


sn 

5(p 



<p]>> 



(A6) 



yields 



pin „> 
> U 12 



0(Gl), 



pin <^~^pin ^> 
U 12 



All neglected terms arising from higher cumulants of 
higher orders of Taylor expansion involve at least two 
propagators G>. Since G> is bounded, terms of order 
G> result in functionals of order d£ 2 for small dl and do 
not contribute to the differential RG. Plugging the con- 
tributions of the last equations into (AC) result in the 
generator (p2|). For further details see also Rcf. |2j3. 



APPENDIX B: KERNELS, COEFFICIENTS 

Here we give explicit expressions for the kernel func- 
tions and the coefficients entering the flow equations. 
This is interesting for illustrative purposes but also im- 
portant to demonstrate, that they are well defined. We 
start from the Gaussian cutoff function ( |l9| ) an d p erform 
calculations in d = 4. We find from Eqs. (pl|) , (38), 
(Egg), and (ph 



5(k) 
a(k) 
c(k) 

P(k) 



3 -fc 2 /2A 2 



1 

5o 
2 

Co 



a(k/v^) 



4A 



2 r 



forfc-0, 
Co In Tj-r- for fc — » 00 . 



-t 2 fc 2 /4A 2 



with go := ,g(x 



0) 



and cq 



4tt 2 7 " u — (47^7) 2 • TneSe 

kernels are apparently well-defined in d = 4 and analytic 
for k 0. 
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FIG. 1. Diagrammatic representation of the flow equations 
(p8[). Black points stand for fields a (x). A shaded oval "ver- 
tex" connection between two points means that they are eval- 
uated at the same position (however, in general their replica 
index is different) and that their difference is the argument 
of a function. A line connecting two fields represents a ker- 
nel function (a broken line represents g; full lines represent 
a or p). Fields connected by a line are confined to the same 
replica. Scaling terms and topologically equivalent diagrams 
are not included in this illustration. At the fixed point an 
oval connection is related to a derivative i?^' . The number of 
lines leaving a vertex equals the order of the derivative. The 
point order is given by the number of vertices. The replica 
order is given by the number of disconnected groups of fields. 
Diagrams in which a g line would connect points that are al- 
ready connected via other lines do not contribute to the flow 
equations since they are of higher order in temperature. If 
both fields of a vertex are linked through lines and an odd 
number of lines leaves the vertex then the diagram vanishes 
for symmetry reasons. 




FIG. 2. Illustration of the fixed-point function R^^cj)) and 
its first four derivatives. Note that d 3 R (1) (0 + ) > 0. The sin- 
gular contribution to d i R^ 1 \(/)) is represented by the vertical 
part of the line near <f> = 0. 
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